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MULTIFRACTAL ANALYSIS BASED ON p-EXPONENTS AND 
LACUNARITY EXPONENTS 

PATRICE ABRY, STEPHANE JAFFARD, ROBERTO LEONARDUZZI, 
CLOTHILDE MELOT, AND HERWIG WENDT 


Abstract. Many examples of signals and images cannot be modeled by lo¬ 
cally bounded functions, so that the standard multifractal analysis, based on 
the Holder exponent, is not feasible. We present a multifractal analysis based 
on another quantity, the p-exponent, which can take arbitrarily large negative 
values. We investigate some mathematical properties of this exponent, and 
show how it allows us to model the idea of “lacunarity” of a singularity at 
a point. We finally adapt the wavelet based multifractal analysis in this set¬ 
ting, and we give applications to a simple mathematical model of multifractal 
processes: Lacunary wavelet series. 


Keywords: Scale Invariance, Fractal, Multifractal, Hausdorff dimension. Hol¬ 
der regularity. Wavelet, Lacunarity exponent, p-exponent 

1. Introduction 

The origin of fractal geometry can be traced back to the quest for non-smooth 
functions, rising from a key question that motivated a large part of the progresses 
in analysis during the nineteenth century: Does a continuous function necessarily 
have points of differentiability? A negative answer to this question was supplied 
by Weierstrass when he built his famous counterexamples, now referred to as the 
Weierstrass functions 


-\-oo 

(1) Wa,b{x) = Y Q”cos(b"7ra:) 

n—0 

where 0 < a < 1, 5 was an odd integer and a 6 > 1 -I- 37r/2. The fact that they are 
continuous and nowhere differentiable was later sharpened by Hardy in a way which 
requires the notion of pointwise Holder regularity, which is the most commonly used 
notion of pointwise regularity in the function setting. We assume in the following 
that the functions or distributions we consider are defined on K. However, most 
results that we will investigate extend to several variables. 

Definition 1. Let f : M —>■ M 6e a loeally bounded funetion, xq € M and let 7 > 0; 
/ belongs to C'^(xo) if there exist C > 0, R > 0 and a polynomial P of degree less 
than 7 such that: 

(2) for a.e. x such that \x — xo\ < R, \f{x) — P{x — xo)\<C\x — xop. 
The Holder exponent of f at xg is 

(3) /i/(xo) = sup {7 : / 


1 


is C^{xg)). 
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The Holder exponent of Wa,b is a constant function, which is equal to H — 

— loga/log6 at every point (see e.g. [13] for a simple, wavelet-based proof); since 
-H < 1 we thus recover the fact that 'Wa,b is nowhere differentiable, but the sharper 
notion of Holder exponent allows us to draw a difference between each of the Weier- 
strass functions, and classify them using a regularity parameter that takes values 
in IR+. The graphs of Weierstrass functions supply important examples of fractal 
sets that still motivate research (the determination of their Hausdorff dimensions 
remains partly open, see m)- In applications, such fractal characteristics have 
been used for classification purposes. For instance, an unorthodox use was the dis¬ 
crimination between Jackson Pollock’s original paintings and fakes using the box 
dimension of the graph supplied by the pixel by pixel values of a high resolution 
photograph of the painting, see P5] . 

The status of everywhere irregular functions was, for a long time, only the one 
of academic counter-examples, such as the Weierstrass functions. This situation 
changed when stochastic processes like Brownian motion (whose Holder exponent 
is J/ = 1/2 everywhere) started to play a key role in the modeling of physical phe¬ 
nomena. Nowadays, experimentally acquired signals that are everywhere irregular 
are prevalent in a multitude of applications, so that the classification and modeling 
of such data has become a key problem. However, the use of a single parameter (e.g. 
the box dimension of the graph) is too reductive as a classification tool in many 
situations that are met in applications. This explains the success of multifractal 
analysis, which is a way to associate a whole collection of fractal-based parameters 
to a function. Its purpose is twofold: on the mathematical side, it allows one to de¬ 
termine the size of the sets of points where a function has a given Holder exponent; 
on the signal processing side, it yields new collections of parameters associated 
to the considered signal and which can be used for classification, model selection, 
or for parameter selection inside a parametric setting. The main advances in the 
subject came from a better understanding of the interactions between these two 
motivations, e.g., see |3] and references therein for recent review papers. 

Despite the fact that multifractal analysis has traditionally been based on the 
Holder exponent, it is not the only characterization of pointwise regularity that can 
be used. Therefore, our goal in the present contribution is to analyze alternative 
pointwise exponents and the information they provide. 

In Section[^we review the possible pointwise exponents of functions, and explain 
in which context each can be used. 

In Section we focus on the p-exponent, derive some of its properties, and 
investigate what information it yields concerning the lacunarity of the local behavior 
of the function near a singularity. 

In Sectionj^we recall the derivation of the multifractal formalism and give appli¬ 
cations to a simple model of a random process which displays multifractal behavior: 
Lacunary wavelet series. 

We conclude with remarks on the relationship between the existence of p-exponents 
and the sparsity of the wavelet expansion. 


This paper partly reviews elements on the p-exponent which are scattered in the 
literature, see e.g. m [i (H HU HO]. New material starts with the introduction 
and analysis of the lacunarity exponent in Section 2.3 the analysis of thin chirps 
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in Section |3.5[ and all following sections, except for the brief reminder on the 
multifractal formalism in Section l4.ll 

2. POINTWISE EXPONENTS 

In this section, unless otherwise specified, we assume that / G An 

important remark concerning the definition of pointwise Holder regularity is that if 
(H) holds (even for 7 < 0 ), then / is bounded in any annulus 0 <r< |a: — XqI < R. 
It follows that, if an estimate such as (H holds for all xq, then / will be locally 
bounded, except perhaps at isolated points. For this reason, one usually assumes 
that the considered function / is (everywhere) locally bounded. It follows that Q 
holds for 7 = 0 so that the Holder exponent is always nonnegative. 

2.1. Uniform Holder regularity. An important issue therefore is to determine 
if the regularity assumption / G is satisfied for real life data. This can be 
done in practice by first determining their uniform Holder exponent, which is 
defined as follows. 

Recall that Lipschitz spaces C®(K) are defined for 0 < s < 1 by 

/ G and 3C, yx,y, \f{x)-f{y)\<C\x-y\y 

If s > 1, they are then defined by recursion on [s] by the condition: / G (^^(M) 
if / G L°° and if its derivative /' (taken in the sense of distributions) belongs 
to C‘^-^(R). If s < 0 , then the C® spaces are composed of distributions, also 
defined by recursion on [s] as follows: / G (^^(K) if / is a derivative (in the sense 
of distributions) of a function g G (^^^^(K). We thus obtain a definition of the 
spaces for any s ^ Z (see [H] for s G Z, which we will however not need to 
consider in the following). A distribution / belongs to if fip G for every 
C°° compactly supported function (p. 

Definition 2. The uniform Holder exponent of a tempered distribution f is 

(4) = sup{s : / G cum- 

This definition does not make any a priori assumption on /: The uniform Holder 
exponent is defined for any tempered distribution, and it can be positive or negative. 
More precisely: 

• If i?™™ > 0, then / is a locally bounded function, 

• if Hy™ < 0 , then / is not a locally bounded function. 

In practice, this exponent is determined through the help of the wavelet co¬ 
efficients of /. By definition, an orthonormal wavelet basis is generated by 
a couple of functions ( 1 ^, "0), which, in our case, will either be in the Schwartz 
class, or smooth and compactly supported (in that case, wavelets are assumed to 
be smoother than the regularity exponent of the considered space). The functions 
ip{x—k), A: G Z, together with 2^/^'0(2'la;—fc), j > 0, A: G Z, form an orthonormal 
basis of L^(M). Thus any function / G L^(K) can be written 

fix) = '^(^k m “ ^) + X! X! '0(2'^ 2^ - A;), 

k j>0 

where the wavelet coefficients of / are given by 

(5) J 


(p{t — k)f{t)dt and Cj^k = 2-’ 


— k)f{t)dt. 
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An important remark is that these formulas also hold in many different functional 
settings (such as the Besov or Sobolev spaces of positive or negative regularity), 
provided that the picked wavelets are smooth enough (and that the integrals ([^ 
are understood as duality products). 

Instead of using the indices (j, fc), we will often use dyadic intervals: Let 

(6) Xi=XU,k)) = 

and, accordingly: c\ = Cj^k and ipxix) = ip{2^x — k). Indexing by dyadic intervals 
will be useful in the sequel because the interval A indicates the localization of the 
corresponding wavelet: When the wavelets are compactly supported, then, 3C > 0 
such that when supp{ip) C [—(7/2, C/2], then supp{tp\) C 2CA. 


A 

23 ’ 23 ) 


In practice, can be derived directly from the wavelet coefficients of / 

through a simple regression in a log-log plot; indeed, it follows from the wavelet 
characterization of the spaces C®, see [22], that: 


( 7 ) 


777*“ = liminf 

>- + CDO 


log ^sup I 


b'A' 


log(2-^) 


This estimation procedure has been studied in more detail in |16j . Three examples 
of its numerical application to real-world functions are provided in Figure 

A multifractal analysis based on the Holder exponent can only be performed 
if / is locally bounded. A way to determine if this is the case consists in first 
checking if > 0. This quantity is perfectly well-defined for mathematical 

functions or stochastic processes; e.g. for Brownian motion, iL™" = 1/2, and 
for Gaussian white noise, 77™“ = —1/2. However the situation may seem less 
clear for experimental signals; indeed any data acquisition device yields a finite 
set of locally averaged quantities, and one may argue that such a finite collection 
of data (which, by construction, is bounded) can indeed be modeled by a locally 
bounded function. This argument can only be turned by revisiting the way that 
Q is computed in practice: Estimation is performed through a linear regression in 
log-log coordinates on the range of scales available in the data and 77™“ can 
indeed be found negative for a finite collection of data. At the modeling level, this 
means that a mathematical model which would display the same linear behavior in 
log-log coordinates at all scales would satisfy 77™“ < 0. 

The quantity 77™“ can be found either positive or negative depending on the 
nature of the application. For instance, velocity turbulence data and price time 
series in finance are found to always have 77™“ > 0, while aggregated count Internet 
traffic time series always have 77™“ < 0. For biomedical applications (cf. e.g., fetal 
heart rate variability) as well as for image processing, 77™“ can commonly be found 
either positive or negative (see Figurej^ P[51 [THl[TTl 125] . This raises the problem of 
using other pointwise regularity exponents that would not require the assumption 
that the data are locally bounded. We now introduce such exponents. 


2.2. The p-exponent for p > 1. The introduction of p-exponents is motivated 
by the necessity of introducing regularity exponents that could be defined even 
when 77™“ is found to be negative; TP{xo) regularity, introduced by A. Calderon 
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and A. Zvffmund in |8], has the advantage of only making the assumption that f 
locally belongs to LP(R). 

Definition 3. Let p > 1 and assume that f G Let a G K; the function f 

belongs to TP{xq) if there exists C and a polynomial of degree less than a such 

that, for r small enough, 

/ 1 pxo+r \ i/p 

(8) y \f{x)-P,,g{x)\Pdxj <Cr°‘. 

Note that the Taylor polynomial P^g of / at xq might depend on p. However, 
one can check that only its degree does (because the best possible a that one can 
pick in Q depends on p so that its integer part may vary with p, see m)- Therefore 
we introduce no such dependency in the notation, which will lead to no ambiguity 
afterwards. 

The p-exponent of / at xq is defined as 

(9) /i^(a;o) = sup{a : / e TP(a;o)}. 

The condition that / locally belongs to LP(K) implies that ([^ holds for a = 
— Ijp, so that h^(xo) > —ifp- 


We will consider in the following “archetypical” pointwise singularities, which 
are simple toy-examples of singularities with a specific behavior at a point. They 
will illustrate the new notions we consider and they will also supply benchmarks on 
which we can compute exactly what these new notions allow us to quantify. These 
toy-examples will be a test for the adequacy between these mathematical notions 
and the intuitive behavior that we expect to quantify. The first (and most simple) 
“archetypical” pointwise singularities are the cusp singularities. 

Let a S K — 2N be such that a > — 1. The cusp of order a at 0 is the function 

(10) C„(x) = 1x1“. 


The case a S 2N is excluded because it leads to a C°° function. However, if a = 2n, 
one can pick 

C2„(x) = 

in order to cover this case also. 

If a > 0, then the cusp Ca is locally bounded and its Holder exponent at 0 is 
well-defined and takes the value a. If a > —1/p, then its p-exponent at 0 is well- 
defined and also takes the value a, as in the Holder case. (Condition a > — 1/p 
is necessary and sufficient to ensure that Ca locally belongs to L^.) Examples for 
cusps with several different values of a are plotted in Figure 

If / e in a neighborhood of Xq for a p > 1, let us define the critical 
Lebesgue index of / at xq by 

(11) Poif) = sup{p : / G L]/c(ffi) in a neighborhood of xq}. 


The importance of this exponent comes from the fact that it tells in practice for 
which values of p a p-exponent based multifractal analysis can be performed. There¬ 
fore, its numerical determination is an important prerequisite that should not be 
bypassed in applications. In Section 3.1 we will extend the definition of po(/) to 
situations where / ^ Lj^^ and show how it can be derived from another quantity, 
the wavelet scaling function, which can be effectively computed on real-life data. 
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2.3. The lacunarity exponent. The p-exponent at xq is defined on the interval 
[l,Po(/)] or [l,po(/)); when the p-exponent does not depend on p on this interval, 
we will say that / has a p-invariant singularity at xg- Thus, cusps are p-invariant 
singularities. 

This first example raises the following question: Is the notion of p-exponent only 
relevant as an extension of the Holder exponent to non-locally bounded functions? 
Or can it take different values with p, even for bounded functions? And, if such is the 
case, how can one characterize the additional information thus supplied? In order 
to answer this question, we introduce a second type of archetypical singularities, 
the lacunary singularities, which will show that the p-exponent may be non¬ 
constant. We first need to recall the geometrical notion of accessibility exponent 
which quantifies the lacunarity of a set at a point, see [TS]. We denote by A4(A) 
the Lebesgue measure of a set A. 

Definition 4. Let H C M. A point xg of the boundary of Ll is a-accessible if there 
exist C > 0 and rp > 0 such that Vr < rg, 

(12) M{nnB{xg,r)) <Cr°^+^. 

The supremum of all values of a such that holds is called the accessibility 
exponent of fl at xg. We will denote it by Sxgi^)- 

Note that is always nonnegative. If it is strictly positive, then H is 

lacunary at xg. The accessibility exponent supplies a way to estimate, through a 
log-log plot regression, the “size” of the part of H which is contained in arbitrarily 
small neighborhoods of xg. The following sets illustrate this notion. 

Let w and 7 be such that 0 < 7 < w; the set Uuj^-y is defined as follows. Let 

(13) 7^,., = [2-^2-“^'+2-^^-]; then = \J 1^^^. 

3>0 

Clearly, at the origin, 

(14) £oiU^,j) = - - 1. 

OJ 

We now construct univariate functions : K —K which permit us to better 
understand the conditions under which p-exponents will differ. These functions will 
have a lacunary support in the sense of Definition 

Let Tp be the Haar wavelet: if = l[o,i/ 2 ) ~ l[i/ 2 ,i) and 

9{x) = 4’i‘2x) — ip{2x — 1 ) 

(so that 9 has the same support as tp but its two first moments vanish). 


Definition 5. Let a G K and 7 > I. The lacunary comb is the function 

00 

(15) F“^(x) = ^ 2-“^0 {2'r^x - 2-^)) . 

i=i 


Note that its singularity is at xg = 0. Numerical examples of lacunary combs 
are provided in Figure 

Note that the support of is so that the accessibility exponent at 0 of 
this support is given by (14). The function F^^^ is locally bounded if and only if 
a > 0. Assume that a < 0; then F“ locally belongs to L^ if and only if a > —^fp. 
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When such is the case, a straightforward computation yields that its p-exponent at 
0 is given by 

(16) (a::o) = - + (-- l) -■ 

"•T- w \UJ / p 

In contradistinction with the cusp case, the p-exponent of at 0 is not a constant 
function of p. Let us see how the variations of the mapping p hF^i^xo) are related 
with the lacunarity of the support of /, in the particular case of We note that 
this mapping is an affine function of the variable q = 1/p (which, in this context, 
is a more natural parameter than p) and that the accessibility exponent of the 
support of can be recovered by a derivative of this mapping with respect to 
q. The next question is to determine the value of q at which this derivative should 
be taken. This toy-example is too simple to give a clue since any value of q would 
lead to the same value for the derivative. We want to find if there is a more natural 
one, which would lead to a canonical definition for the lacunarity exponent. It is 
possible to settle this point through the following simple perturbation argument: 
Consider a new singularity F that would be the sum of two functions Fi = F/// 
and F 2 = F///l^^ with 

(17) 0 < «! < a 2 and 71 > 72 . 


The p-exponent of F (now expressed in the q variable, where <7 = 1 /p) is given by 


(18) 


q hf (xq) = min 


OJ 


\ OJ ) U! \ UJ J 


The formula for the lacunarity exponent should yield the lacunarity of the most 
irregular component of F\ since F € the Holder exponent is the natural way 
to measure this irregularity. In this respect, the most irregular compo nen t is Fi; 
the lacunarity exponent should thus take the value (^ — l)- But, since (17) allows 
the shift in slope of the function (18) from — l) to — l) to take place at a 
q arbitrarily close to 0 , the only way to obtain this desired result in any case is to 
pick the derivative of the mapping q —>■ precisely at g = 0 . 


A similar perturbation argument can be developed if po(/) < 00 with the con¬ 
clusion that the derivative should be estimated at the smallest possible value of g, 
i.e. for 

hence the following definition of the lacunarity exponent. 


Definition 6. Let f € in a neighborhood of Xg for a p > 1, and assume 
that the p-exponent of f is finite in a left neighborhood o/po(/). The laeunarity 
exponent of f at xg is 

d 


(19) 


c,M = ^ (t 


1/9 

/ 


9=9o(/) + 


Remarks: 

• Even if the p-exponent is not defined at po(/), nonetheless, because of the 


concavity of the mapping g —>■ h^J‘‘(xg) (see Proposition 
right derivative is always well-defined, possibly as a limit. 


3.1 


below), its 
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• As expected, the lacunarity exponent of a cusp vanishes, whereas the lacu- 
narity exponent of a lacunary comb coincides with the accessibility expo¬ 
nent of its support. 

• The condition £y(a;o) ^ 0 does not mean that the support of / (or of / —P) 
has a positive accessibility exponent (think of the function -I- g where 
g is a C°° but nowhere polynomial function). 

• The definition supplied by (1^ bears similarity with the definition of the 
oscillation exponent (see [31 I16j and ref. therein) which is also defined 
through a derivative of a pointwise exponent; but the variable with respect 
to which the derivative is computed is the order of a fractional integration. 
The relationships between these two exponents will be investigated in a 
forthcoming paper [20] • 


3. Properties of the p-exponent 


In signal and image processing, one often meets data that cannot be modeled 
by functions / S see Figure It is therefore necessary to set the analysis in 
a wider functional setting, and therefore to extend the notion of TP{xq) regularity 
to the case p <1. 

3.1. The case p < 1. The standard way to perform this extension is to consider 
exponents in the setting of the real Hardy spaces (with p < 1) instead of 
spaces, see miin]. First, we need to extend the definitions that we gave to the 
range p € (0,1]. The simplest way is to start with the wavelet characterization of 
LP spaces, which we now recall. 

We denote indifferently by Xj,k or the characteristic function of the interval 
A (= Xj^k) defined by (§. The wavelet square function of / is 

1/2 


\^3X\^X3,k{x) 

.U,k)ez^ 


Then, for p > I, 

( 20 ) 


/ G LP(R) 


{Wf{x))^ dx < 00 , 


see [22]. The quantity (/ {Wf{x)Y dxY^^ is thus equivalent to || / ||p. One can 


then take the characterization supplied by (20) when p > 1 as a definition of the 


Hardy space (when p < 1); note that this definition yields equivalent quantities 
when the (smooth enough) wavelet basis is changed, see |22|. This justifies the 
fact that we will often denote by the space , which will lead to no confusion; 
indeed, when p < 1 this notation will refer to Rp, and, when p > I it will refer to 
LP. 


Note that, if p = I, (20) does not characterize the space but a strict subspace 


of (the real Hardy space which consists of functions of whose Hilbert 
transform also belongs to L^, see [22]). 

Most results proved for the LP setting will extend without modification to the 
RP setting. In particular, TP regularity can be extended to the case p < 1 and has 
the same wavelet characterization, see [21] . All definitions introduced previously 
therefore extend to this setting. 
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The definition of T^{xq) regularity given by Q is a size estimate of an norm 
restricted to intervals [xg — r, Xg + r]. Since the elements of can be distributions, 
the restriction of / to an interval cannot be done directly (multiplying a distribution 
by a non-smooth function, such as a characteristic function, does not always make 
sense). This problem can be solved as follows: If / is an open interval, one defines 
II / \\hp(i) = inf II 9 lip; where the infimum is taken on the g G such that f = g 
on /. The TP condition for p < 1 is then defined by: 

/ e TP(xg) ^ II / ||pp.((,„_,,,„+,))< C 

also when p < 1. We will show below that the p-exponent takes values in [— 1 /p, -|-oo]. 


3.2. When can one use p-exponents? We already mentioned that, in order 
to use the Holder exponent as a way to measure pointwise regularity, we need 
to check that the data are locally bounded, a condition which is implied by the 
criterion > 0, which is therefore used as a practical prerequisite. Similarly, 

in order to use a p-exponent based multifractal analysis, we need to check that 
the data locally belong to LP or HP, a condition which can be verified in practice 
through the computation of the wavelet scaling function, which we now recall. 

The Sobolev space is defined by 

Vs G K, Vp > 0, / G LP’" ^ {Id - A)"/V g LP, 

where the operator {Id — is the Fourier multiplier by (1 -I- ^md we 

recall our convention that LP denotes the space HP when p < 1, so that Sobolev 
spaces are defined also for p < 1. 


Definition 7. Let f be a tempered distribution. The wavelet scaling function of f 
is defined by 

(21) Vp > 0, r]f{p) = p sup{s : / G LP’®}. 


Thus, Vp > 0: 

• If Vfip) > 0 then / G Lf^^. 

• If Vfip) < 0 then / Lf^^. 

The wavelet characterization of Sobolev spaces implies that the wavelet scaling 
function can be expressed as (cf. [TT] 1 


( 22 ) 


Vp > 0, 


pj (p) = lim inf 

j—>-+oo 


log 



log(2-J) 


This provides a practical criterion for determining if data locally belong to IP, 
supplied by the condition rif{p) > 0. The following bounds for pg(/) follow: 


sup{p : r]f{p) > 0} < po{f) < inf{p : p/(p) < 0}, 


which (except in the very particular cases where pf vanishes identically on an 
interval) yields the exact value of pg(/). 

In applications, data with very different values of pg(/) show up; therefore, in 
practice, the mathematical framework supplied by the whole range of p is relevant. 
As an illustration, three examples of real-world images with positive and negative 
uniform Holder exponents and with critical Lebesgue indices above and below pg = 
1 are analyzed in Figure 
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Figure 1. Real-world images (top row) of Romanesco broc¬ 
coli (left column), fern leaves (center column) and a patch of a 
hyperspectral image of the Moffett field, acquired by the AVIRIS 
instrument (spectral band 90, right column). Bottom row: corre¬ 
sponding numerical estimation of uniform Holder exponents R™*”, 
wavelet scaling functions r]f{p) and critical Lebesgue indices po, re¬ 
spectively. 


3.3. Wavelet characterization of p-exponents. In order to compute and prove 
properties of p-exponents we will need the exact wavelet characterization of 
see [min]. Let A be a dyadic interval; 3A will denote the interval of same center 
and three times wider (it is the union of A and its two closest neighbors). For 
Xq € denote by Xj{xo) the dyadic cube of width 2“^ which contains Xq- The 
local square functions at a;o are the sequences defined for j > 0 by 

1/2 

I CApXA(a;) 

1 AC3Aj(a:o) 




Recall that (cf. [2T] 1 

(23) / G TP{xo) if and only if 3C > 0, Vj > 0 


yvi 




The following result is required for the definition of the lacunarity exponent in 
(191 to make sense, and implies that Definition]^ also makes sense whenpo(/) < 1- 

Proposition 3.1. Let p,q & (0,-l-oo], and suppose that f G TP{xo) H T^{xo); let 
9 G [0,1]. Then f G T^{xo), where 

- = —I - and 'y = 9a + il — 9)/3. 


It follows that the mapping q —>■ h^j‘^{xif) is concave on its domain of definition. 


Proof: When p,q < oo, the result is a consequence of (23). Holder’s inequality 
implies that 





e/p 


Wf 

f,Xo 

< 

r 

wi 

f,^o 

p 

Wi 

f,xo 


g 
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Figure 2. Cusps with exponents a = {+0.3, —0.2, —2} (from 
top to bottom row, respectively): functions (left column) and esti¬ 
mation of p-exponents and lacunarity exponents (center and right 
column, respectively). The critical Lebesgue indices are given by 
Po = 1+00, 5, 0.5}, respectively. 


We thus obtain the result for p, g < oo. The case when p or g = +oo does not 
follow, because there exists no exact wavelet characterization of C°‘(xo) = T^(xo); 
however, when p,g > 1, one can use the initial definition of TP(xq) and C'“(a;o) 
through local and L°° norms and the result also follows from Holder’s inequal¬ 
ity; hence Proposition |3.1| holds. 


If / G then || Wf ||p< C. Since < W/, it follows that 




lp< c, 


so that (231 holds with a = —1/p. Thus p-exponents are always larger than — 1/p 
(which extends to the range p < 1 the result already mentioned for p > 1). Note 
that this bound is compatible with the existence of singularities of arbitrary large 
negative order (by picking p close to 0). The example of cusps will now show that 
the p-exponent can indeed take values down to — 1/p. 

3.4. Computation of p-exponents for cusps. Typical examples of distributions 


for which the p-exponent is constant (see Proposition 3.2 below) and equal to a given 
value a < — 1 are supplied by the cusps Ca, whose definition can be extended to 
the range a < —1 as follows: First, note that cusps cannot be defined directly 
for a < —1 by (10) because they do not belong to so that they would be ill- 
defined even in the setting of distributions (their integral against a compactly 
supported function p may diverge). Instead, we use the fact that, if a > 1, then 
C" = a(a — l)Ca- 2 ! which indicates a way to define by recursion the cusps Ca, when 
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Figure 3. Lacunary combs with po = +(X) (top row) and po = 
9.1 (bottom row): functions (left column) and estimation of p- 
exponents and lacunarity exponents (center and right column, re¬ 
spectively). 


a < — 1 and a as follows: 

if a < 0, 


C — 

^OL - 


1 


-C' 


(a-Hl)(a-h2)''“+2’ 

where the derivative is taken in the sense of distributions. The Ca are thus defined 
as distributions when a is not a negative integer. It can also be done when a is a 
negative integer, using the following definition for a = 0 and —1: 


Co = log(|a:|) and 
where P.V. stands for “principal value”. 


C_i =Co = P.V. 


Proposition 3.2. If a > 0, the cusp Ca belongs to and its p-exponent is a. If 
a <0, the cusp Ca belongs to for p < —1/a and its p-exponent is a. 

Proof of Proposition |3.2t The case a > 0 and p > 1 has already been consid¬ 
ered in [ini . In this case, the computation of the p-exponent is straightforward. 

Note that, when a G (—1, 0) and p > 1 the computations are similar. We thus focus 
on the distribution case, i.e. when p < 1. The global and pointwise regularity will 
be determined through an estimation of the wavelet coefficients of the cusp. We 
use a smooth enough, compactly supported wavelet basis and we denote by c^g^ the 
wavelet coefficients of the cusp 

^j.k — {f^j.klCa} ■ 

The selfsimilarity of the cusp implies that 

(24) Vj, fc cyfe = 2““-^co,fe; 

additionally, as soon as k is large enough so that the support of ip{x — k) does not 
intersect the origin, the cusp is (7°° in the support of ip{x — k) and coincides with 
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the function |a;|“. An integration by parts then yields that, for any N smaller than 
the global regularity of the wavelet, 


co,fc = [— 
so that the sequence co,fc satisfies 


(— 1)^ J — k) a{a — 1) • • • (a — N)\x 


a-N 


dx, 


(25) 


|co,fc| < 


Cat 


{l + \k\r 

where N can be picked arbitrarily large. The estimation of the norm of the 


wavelet square function follows easily from (24) and (25), and so does the lower 


bound for the p-exponent. The upper bound is obtained by noticing that one of 
the co,fc does not vanish (otherwise, all Cj^k would vanish, and the cusp would be a 
smooth function at the origin). Therefore, there exists at least one fcg such that Vj, 
Cj^ko = C2~°‘^, and the wavelet characterization of regularity then yields that 
hP{xo) < a. 

Three examples of cusps and numerical estimates of their p-exponents and lacu- 
narity exponents are plotted in Figure 

3.5. Wavelet characterization and thin chirps. In practice, we will derive 
regularity from simpler quantities than the local square functions. The p-leaders 
of / are defined by local P norms of wavelet coefficients as follows: 

i/p 

(26) ' ~ 


\A'C3A / 


(they are finite if / € L(),^(M‘^), see [IS])- Note that, if p = +oo, the corresponding 
quantity is usually denoted by d\ and simply called the wavelet leaders; we 
have 


(27) 


dx := = sup |cA'|- 

A'C3A 


The notion of TP regularity can be related to p-leader coefficients (see [miia 

my- 


If T]f{p) > 0, then h^Jxo) = liminf 

-} —^ -I- oo 




^ ^ '/log(2 

Our purpose in this section is to introduce new “archetypical” pointwise singu¬ 
larities which will yield examples where the p-exponent and the lacunarity exponent 


can take arbitrary values. Because of (28), it is easier to work with examples that 
are defined directly by their wavelet coefficients on a smooth wavelet basis. We 
therefore develop new examples rather than extending the lacunary combs of Sec¬ 
tion [2]3| 

Definition 8. Let a,b G (0,1) satisfying 0 < b < 1 — a, and let a S K. The thin 
chirp Ta,b,a is defined by its wavelet series 

f~a^b,OL ^ ^ 

j>0 fcGZ 
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Figure 4. Thin chirps withpo = oo (top row) andpo = 3.2 (bot¬ 
tom row): functions (left column) and estimation of p-exponents 
and lacunarity exponents (center and right column, respectively). 


where 

=2-“J if fc G [2(1-“)^, 2(1-“)^ + 2^^] 
= 0 otherwise. 


The following results are straightforward, using the wavelet characterization of 
LP and regularity. 


Proposition 3.3. The thin chirp Ta,b,a is bounded if and only if a > 0. 

If a < 0, Po(Ta,b,a) = —-■ 


The p-exponent ofTa,b,a at the origin is 





a 

a 


Note that, if the wavelets are compactly supported, then for j large enough the 
pack of 2^-1 successive wavelets with non-vanishing coefficients covers an interval of 
length 2“-'2^-1 at a distance 2““^ from the origin, so that the accessibility exponent 
of the support of Ta,b,a is (1 — a — &)/a: Thus, it coincides with the lacunarity 
exponent of Ta,b,a as expected. 


Illustrations of thin chirps and the numerical estimation of their p-exponents and 
lacunarity exponents are provided in Figure 

3.6. p-exponent analysis of measures. Several types of measures (such as mul¬ 
tiplicative cascades) played a central role in the development of multifractal anal¬ 
ysis. Since measures (usually) are not functions, their p-exponent for p > 1 is 
not defined. Therefore, it is natural to wonder if it can be the case when p < 1. 
This is one of the purposes of Proposition |3.4[ which yields sufficient conditions 
under which a measure p satisfies rj^ip) > 0 for p < 1, which will imply that its 
p-exponent multifractal analysis can be performed. An important by-product of 
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using p-exponents for p < 1 is that it offers a common setting to treat pointwise 
regularity of measures and functions. 

Recall that dimB(^) denotes the upper box dimension of the set A. 


Proposition 3.4. Let ^ be a measure; then its wavelet scaling function satisfies 
Furthermore, if fj, does not have a density whieh is an function, then 

Vui^) = 0- 

Additionally, if g, is a singular measure whose support supp(p) satisfies 


(29) 


:= dimB{supp{p)) < 1, 

then 



(30) 

Vp < 1, 

VuiP) > (1-<5m)(1-p) 

and 



(31) 

Vp > 1, 

VuiP) < (1-<^m)(1-p) 


Remarks: 


(30) expresses the fact that, if p, has a small support, then its Sobolev 
regularity is increased for p < 1. This is somehow counterintuitive, since 
one expects a measure to become more singular when the size of its support 


shrinks; on the other hand (31) expresses that this is actually the case when 

p > 1. 

Condition < 1 is satisfied if p is supported by a Cantor-like set, or by a 
selfsimilar set satisfying Hutchinson’s open set condition. 


(30) has an important consequence for the multifractal analysis of measures: 


Indeed, if < 1, then p^(p) > 0 for p < 1, so that the classical mathemat¬ 
ical results concerning the multifractal analysis based on the p-exponent 
apply, see Section]^ 

A slightly different problem was addressed by H. Triebel: In [57], he de¬ 
termined under which conditions the scaling functions commonly used in 


the multifractal analysis of probability measures (see (37) below) can be 


recovered through Besov or Triebel-Lizorkin norms (or semi-norms). 
Proof of Proposition [374] If p is a measure, then for any continuous bounded 
function / 


(32) \{p\f)\ < C \\ f \U . 

We pick 

f = where = ±1, 

k 

so that / is continuous and satisfies || / ||oo< C, where C depends only on the 
wavelet (but not on the choice of the £j,fc). Denoting by Cj^k the wavelet coefficients 
of p, we have 

(mI/) — ^ ( ^j,k I '4^j,kdp — 2 ^ ( ^j,k^j,k- 

k k 

Picking Ej^k = sgn(cj_/c) it follows from ( [^ that 

(33) 2-^'^|c,-fc|<C, 
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or, in other words, fj, belongs to the Besov space , which implies that ?7^(1) > 0, 

see [131 [21] ■ 

On other hand, if /i ^ then using the interpretation of the scaling function 
in terms of Sobolev spaces given by (21), we obtain that ?7^(1) < 0. Hence the first 
part of the proposition holds. 

We now prove (30). We assume that the used wavelet is compactly supported, 
and that its support is included in the interval [—2^,2*] for an Z > 0 (we pick the 
smallest I such that this is possible). Let S > dimB{supp{fi))-, for j large enough, 
supp{p) is included in at most 2[^^1 intervals of length 2“-^ . It follows that, at scale 
j, there exist at most •2-2^ wavelets (tpj,k)kez whose support intersects the 
support of p. Thus for j large enough, there are at most wavelet coefficients 
that do not vanish. 

Let p e (0,1), q = 1/p and r be the conjugate exponent of q, i.e. such that 
1/g + 1/r = 1. Using Holder’s inequality, 

k \ k / \ k / 

where the sums are over at most C2^^ terms; thus 

Eic..fcr< 


Using (33), we obtain that 






so that r]fi{p) > (1 — (5)(1 — p)- Since this is true Vd > <5^, (30) follows. 

We now prove (31). Let p > 1 and let q be the conjugate exponent. Using 
Holder’s inequality, 


E icj.fei ^ E 


i/p 


El' 


1/9 


Let again S > Sfj,; using the fact that the sums bear on at most 2^^ terms, and that 
the left-hand side is larger than (72/, we obtain that 

>C2n-^^i\ 


which can be rewritten 


2 -^ E ^ ^ 2-^2P^2-'^/p/i, 


so that r}^{p) < (1 — p){l — 6)] since this is true V5 > 5^, (31) follows, and Propo¬ 
sition |3.4| is completely proved. 


Since p = 1 is a borderline case for the use of the 1-exponent one may expect that 
picking p < 1 would yield r]^{p) > 0 (in which case one would be on the safe side 
in order to recover mathematical results concerning the p-spectrum, see [HIT]). 
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However, this is not the case, since there exist even continuous functions / that 
satisfy yp > 0, r]f(p) = 0. An example is supplied by 

j>0 fcGZ 

4. Multifractal analysis of lacunary wavelet series 


Multifractal analysis is motivated by the observation that many mathematical 
models have an extremely erratic pointwise regularity exponent which jumps every¬ 
where; this is the case e.g. of multiplicative cascades or of Levy processes, whose 
exponents h satisfy that 


(34) 


a.s. Vxo, limsup/i(a;) — liminf h(a;) 

X^Xq X^Xq 


is bounded from below by a fixed positive quantity (we will see that this is also the 
case for lacunary wavelet series). This clearly excludes the possibility of any robust 
direct estimations of h. The driving idea of multifractal analysis is that one should 
rather focus on alternative quantities that 


• are numerically computable on real life data in a stable way, 

• yield information on the erratic behavior of the pointwise exponent. 

Furthermore, for standard random models (such as the ones mentioned above) we 
require these quantities not to be random (i.e. not to depend on the sample path 
which is observed) but to depend on the characteristic parameters of the model 
only. The relationship between the multifractal spectrum and scaling functions 
(initially pointed out by U. Frisch and G. Parisi in [531; below) satisfies 

these requirements. 

We now recall the notion of multifractal spectrum. We denote by dim(A) the 
Hausdorff dimension of the set A. 


Definition 9. Let h{x) denote a pointwise exponent. The multifractal spectrum 
d{H) associated with this pointwise exponent is 

d{H) = dim{a; : h{x) = H}. 


In the case of the p-exponent, the sets of points with a given p-exponent will be 
denoted by F^{H): 

(35) FPiH) = {xo: hPixo)=H}, 

and the corresponding multifractal spectrum (referred to as the p-spectrum) is 
denoted by d^{H)-, in the case of the lacunarity exponent, we denote it by d^{L). 


4.1. Derivation of the multifractal formalism. We now recall how d{F[) is 
expected to be recovered from global quantities effectively computable on real-life 
signals (following the seminal work of G. Parisi and U. Frisch [23] and its wavelet 
leader reinterpetation m)- A key assumption is that this exponent can be derived 
from nonnegative quantities (which we denote either by ej^k or e;^), which are 
defined on the set of dyadic intervals, by a log-log plot regression: 


(36) 


h(xo) = liminf 


log jexpxo)) 

log(2-f) 


It is for instance the case of the p-exponent, as stated in (23) or (28), for which the 
quantities e\ are given by the p-leaders d^. 
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In the case of the lacunarity exponent, quantities e\ can be derived as follows: 
Let Ag > 0 small enough be given. If / has a I/g-exponent H and a lacunarity 
exponent L at Xq then its 1/g-leaders satisfy 

dy’^(xo) ~ 2 -^^ 

and its l/(q + Ag)-leaders satisfy 

^i/(9+A9)(^o) ^ 

we can eliminate from these two quantities by considering the /1-leaders: 


,l/(g+Ag)\ 1/Ag 


4 ■■ = 




l-Lj 


(this argument follows a similar one developed in |16l Ch. 4.3] for the derivation of 
a multifractal analysis associated with the oscillation exponent). 


The multifractal spectrum will be derived from the following quantities, referred 
to as the structure functions, which are similar to the ones that come up in the 


characterization of the wavelet scaling function in (22): 


SfirJ) = {‘2 




The scaling function associated with the collection of {e\) is 

log(S'/(r,j)) 


(37) 


Vr e 


ff{r) = liminf 

ji—>-+oo 


log(2-fy 


Let us now sketch the heuristic derivation of the multifractal formalism; (37) means 
that, for large j, 

5/(r,j)~2-««L 

Let us estimate the contribution to Sf{r,j) of the dyadic intervals A that cover the 
points of Eh- By definition of Eh, they satisfy e\ ~ 2“'^-^; by definition of d{H), 
since we use cubes of the same width 2~^ to cover Eh, we need about such 

cubes; therefore the corresponding contribution is of the order of magnitude of 


When j 
that 

(38) 


+ 00 , the dominant contribution comes from the smallest exponent, so 
((r) = inf(l — d(H) + Hr). 

H 

By construction, the scaling function f[r) is a concave function on M, see 


[24] which is in agreement with the fact that the right-hand side of (38) necessarily 
is a concave function (as an infimum of a family of linear functions) no matter 
whether d{H) is concave or not. If d{H) also is a concave function, then the 


Legendre transform in (38) can be inverted (as a consequence of the duality of 


convex functions), which justifies the following assertion. 

Definition 10. A nonnegative sequence {e\), defined on the dyadic intervals, fol¬ 
lows the multifractal formalism if the associated multifractal spectrum d{H) satisfies 


( 39 ) 


d(H) = inf (1 — C(r) -f Hr). 

rGR 








MULTIFRACTAL ANALYSIS BASED ON p-EXPONENTS AND LACUNARITY EXPONENTSI9 


The derivation given above is not a mathematical proof, and the determination 


of the range of validity of (39) (and of its variants) is one of the main mathematical 


problems concerning multifractal analysis. If it does not hold in complete generality, 
the multifractal formalism nevertheless yields an upper bound of the spectrum of 
singularities, see [23l [131 [24] : As soon as (36) holds. 


d(H) < inf (1 — ((r) + Hr), 
res. 

In applications, multifractal analysis is often used only as a classification tool 
in order to discriminate between several types of signals; then, one is not directly 


concerned with the validity of (39) but only with a precise computation of the 
new multifractal parameters supplied by the scaling function, or equivalently its 
Legendre transform. Note that studies of multifractality for the p-exponent have 
been performed by A. Fraysse who proved genericity results of multifractality for 
functions in Besov or Sobolev spaces in m- 


4.2. Description of the model and global regularity. In this section, we ex¬ 
tend to possibly negative exponents the model of lacunary wavelet series introduced 
in [12]. We assume that ^ is a wavelet in the Schwartz class (see however the re¬ 
mark after Theorem [^ which gives sufficient conditions of validity of the results of 
this section when wavelets of limited regularity are used). Lacunary wavelet series 
depend on a lacunarity parameter j] G (0,1) and a regularity parameter 
a G K. At each scale j > 0, the process Xa.y has exactly [2'^^] nonvanishing wavelet 
coefficients on each interval [/, Z -I- 1) (Z G Z), their common size is 2~°‘^, and their 
locations are picked at random: In each interval [Z,Z -|- 1) (Z G Z), all drawings of 
[2"^^] among the 2^ possibilities ^ G [Z,Z -|- 1) have the same probability. Such a 
series is called a lacunary wavelet series of parameters (a, rj). Note that, since a 
can be arbitrarily negative, Xa^n can actually be a random distribution of arbitrary 
large order. By construction 


jjrnin 


= a, 


and, more precisely, the sample paths of Xa.n are locally bounded if and only if 
a > 0. The case considered in [T2| dealt with a > 0, and was restricted to the 
computation of Holder exponents. Considering p-exponents allows us to extend the 
model to negative values of a, and also to see how the global sparsity of the wavelet 
expansion (most wavelet coefficients vanish) is related with the pointwise lacunarity 
of the sample paths. Note that extensions of this model in different directions have 
been worked out in 

Since we are interested in local properties of the process X, we restrict our 
analysis to the interval [0,1) (the results proved in the following clearly do not 
depend on the particular interval which is picked); we can therefore assume that 
fcG {0,---2J-1}. 

We first determine how a and r\ are related with the global regularity of the 


sample paths. The characterization (22) implies that the wavelet scaling function 
is given by 


(40) 

It follows that 


Vp > 0, ?7jc„,,(p) = ap-?? + 1. 


Po ■= Po{Xa,rt) 



if a < 0 
if a > 0. 
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Note that po always exists and is strictly positive, even if a takes arbitrarily large 
negative values. We recover the fact that p-exponents allow us to deal with sin¬ 
gularities of arbitrarily large negative order. We will see that this is a particular 
occurrence of a general result, see Proposition |5.1| the key property here is the 
sparsity of the wavelet series. 

4.3. Estimation of the p-leaders of Xa.,n- An important step in the determi¬ 
nation of the p-exponent of sample paths of ^ at every point is the estimation 
of their p-leaders. We now assume that p < po, so that the sample paths of Xa^n 
locally belong to and the p-exponent of Xa^n is well-defined everywhere. Recall 
that the p-leaders are defined by 



The derivation of the p-exponent of everywhere will be deduced from the 
estimation of the size of the p-leaders of Xa^^^- A key result is supplied by the 
following proposition, which states that the size of the p-leaders of a lacunary 
wavelet series is correctly estimated by the size of the first nonvanishing wavelet 
coefficient of smaller scale that is met in the set {A' : A' C 3A}. 

Proposition 4.1. Let a S M, p S (0,1) and let X^^ri be a laeunary wavelet series of 
parameters (a, p); for each dyadic interval X (of width 2~^), we define j' ( = j'{X)) 
as the smallest random integer such that 

3A^ C 3A such that |A^| = 2~^ and cy 0. 

Then, a.s. 3J, 3C, C' > 0 such that Vj > J, VA of scale j 

<l^< 

Proof: This result will be implied by the exponential decay rate 2~^^ that 
appears in the definition of p-leaders together with the lacunarity of the construc¬ 
tion; we will show that exceptional situations where this would not be true (as a 
consequence of local accumulations of nonvanishing coefficients) have a small prob¬ 
ability and ultimately will be excluded by a Borel-Cantelli type argument. We now 
make this argument precise. For that purpose, we will need to show that the spar¬ 
sity of wavelet coefficients is uniform, which will be expressed by a uniform estimate 
on the maximal number of nonvanishing coefficients cy that can be found for (at 
a given scale j') included in a given interval 3A. Such an estimate can be derived 
by interpreting the choice of the nonvanishing wavelet coefficients in the construc¬ 
tion of the model as a coarsening (on the dyadic grid) of an empirical process. 
Let us now recall this notion, and the standard estimate on the increments of the 
empirical process that we will need. 

Let Nj = [2’^'f] denote the number of nonvanishing wavelet coefficients at scale 
j. We can consider that the corresponding dyadic intervals A have been obtained 
first by picking at random Nj points in the interval [0,1] (these points are now 
Nj independent uniformly distributed random variables on [0,1]), and then by 
associating to each point the unique dyadic interval of scale j to which it belongs. 
Let Pf be the process starting from 0 at t = 0, which is piecewise constant and 
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which jumps by 1 at each random point thus determined. The family of processes 

(42) = 

is called an empirical process on [0,1]. The size of the increments of the empirical 
process on a given interval yields information on the number of random points 
picked in this interval. If it is of length Z, then the expected number of points is 
l[2^^], and the deviation from this average can be uniformly bounded using the 
following result of W. Stute which is a particular case of Lemma 2.4 of |26) . 

Lemma 4.2. There exist two positive constants C[ and such that, ifO < I < 1/8, 
Njl > 1 and 8 < A < C[ \jNjl, 

P I sup \a{ — ai\ > A-\/Z I < ^e~^ . 

\\t-s\<l J « 


Rewritten in terms of P/, this means that 


(43) P ( sup |P/ - Pi - N,{t - s)| > A^] < 

\\t-s\<l J i 

Recall that the assumption A' C 3A implies that 3 • 2~^ > 2~^ . We will apply 


Lemma 4.2 differently for small values of j' where the expected number of nonva¬ 
nishing coefficients ca' that can be found for A' (at a given scale /) included in a 
given interval A is very small, and the case of large j' where this number increases 
geometrically. 


2"^' > 7^2"^/''. 
/2r\ —r 


We first assume that 

( 44 ) 2 -^ > r 

We pick intervals of length I = j' 2~'^^ and, for the constant A in State’s lemma, 
we pick A = j. Then (43) applied with N = [2'^^ ] yields that, with probability at 
least 1 — e~^ , the number of intervals A' of scale j' picked in such intervals is 

2^U + Oif)=Oif). 

We now assume that 


(45) 


2-^' <f 2-^/T 


Then we pick intervals of length 1 = 3-2 ^, and A = j + j'. Then (43) applied 
with N = [2^^ ] yields that, with probability at least 1 — ^ , the number of 

intervals A' of scale j' picked in such intervals is 

(46) 2^^'I + 0{{j + < 2 ■ 2’'^'Z. 

We are now ready to estimate the size of l\, assuming that all events described 
above happen (indeed, we note that the probabilities such that these events do not 
happen have a finite sum, so that, by the Borel-Cantelli lemma, they a.s. all occur 
for j large enough). 

At scales j' which satisfy (44), if at least one of the A' C 3A does not vanish, 
th en t here are at most of them, and the corresponding contribution to the sum 
in (41) lies between \{c\i)p2~^^~^ ^ | and /^|(ca')*’ 2“*-'^“'^ ^ |. At scales j' which satisfy 
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( [ 45 ] ) , the contribution of the wavelet coefficients of scale f to the sum lies between 
and its double. Since cy = , the condition p < po implies 

that these quantities decay geometrically, so that the order of magnitude of the 


p-leader is given by the first non-vanishing term in the sum. Hence Proposition 4.1 
holds. 


4.4. p-exponents and lacunarity. We now derive the consequences of Proposi¬ 
tion 4.1 for the determination of the p-exponents of Xa^rj at every point. We first 


determine the range of p-exponents. First, note that all p-leaders have size at most 
2““-^, so that the p-exponent is everywhere larger than a. In the opposite direction, 
as a consequence of (461, every interval 3A of scale j includes at least one nonvan¬ 
ishing wavelet coefficient at scale j/rj -I- (log j)^; therefore, all p-leaders have size at 
least 

2-a(i-|-log j)^-i(i-i + (logi)^)^ 

It follows that the p-exponents are everywhere smaller than 


(47) 

We have thus obtained that 

a.s. Vp < Po, Vxo S 


a 


a<h\ (a;o) < i7„ 


For each j, let denote the subset of [0,1] composed of intervals 3A (A S Aj) 
inside which the first nonvanishing wavelet coefficient is attained at a scale I < [wj], 
and let 

= lim sup . 

Proposition |4.1| implies that, if Xq ^ E,^, then, for j large enough, all wavelet leaders 
hj{xo) are bounded by 

■2rt — Oi^— - ( 

J Z ^ P \ V ^ 

SO that: 

(48) if xq 4 Ei^, then (xq) > auj + — - 

p 

On other hand, if xq € E^, then there exists an infinite number of p-leaders l\ (xg) 
larger than 

cy — a.^ — - ( ^ — 

Z V p\v ^ 

SO that: 

( 49 ) if xo G Ei^, then (xq) < acu H-. 

p 


It follows from (48) and (49) that the sets of points where the p-exponent takes the 
value 

uj — 1 

U = auj -\ - 

P 

are the sets 

H^= f] E^,- \J E^,. 

UJ'^UJ Uj'<.(jJ 

We have thus obtained the following result. 
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Proposition 4.3. Let a G K, 17 € (0,1) and Let Xa^n be a lacunary wavelet series 
of parameters (a, rj). Let p < Pq; the sets of points with a given p-exponent are the 
sets 

= for - = 

and additionally, if xq S Hi^, then 

^x^,^{xo) = w - 1. 

Remark: We actually do not need the wavelet used to be in the Schwartz 
class for Theorem to be true. One can verify that, if the uniform regularity of 
the wavelet is larger than max(|a|, |Rmax|)) then all previous computations remain 
valid. 


In order to determine the p-spectra and the lacunarity spectrum, one has to 
determine the Hausdorff dimensions of the sets . We note that these sets do not 
depend on a and on p, but only on the parameter iw and on the random drawing of 
the locations of the non-vanishing wavelet coefficients. When a > 0, the dimensions 
of these sets (expressed in a slightly different way) were determined in [T^], where 
it is shown that 

dim(Rij) = pui. 

The following result follows. 


Theorem 1. Let a G M, ly S (0,1) and let Xa,n be a lacunary wavelet series of 
parameters ( 0 , 77 ); the p-speetrum of Xa,n is supported by the interval [a,Hmax\ 
and, on this interval, 

a.s. yp<po, VH, dP{H) = ■ 

a + 1/p 

Furthermore, its lacunarity spectrum is given by 

a.s. VL € [0, 1/77 — 1], d^{L) = ^{L + 1). 


Remark: It is also shown in [12] that all the sets are everywhere dense, so 
that the quantity (34) is equal everywhere to H^ax — O'. 


For the sake of completeness, we now sketch how these dimensions can be com¬ 
puted. We start by estimating the size of Note that the number of intervals 
3A which comprise E/, is bounded by 

-p H-h < C'2''“T 

Using these intervals for j > J as an e-covering, we obtain the following bound for 
the Hausdorff dimension of E^ 


(50) 


dim(i?,^) < rjUJ. 


We now consider the sets it follows from 


and (49) that 

= n - u 


Since Vw' < u, F[^ C E^/, it follows from (49) that 


dim(Rij) < rjUj. 
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In order to get a lower bound on the Hausdorff dimension of we will need 
the following (slightly) modified notion of 5-dimensional Hausdorff measure. 

Definition 11. Let H C M. For e > 0 and S S [0,1], let 

M,^’^(H)=inf (^^|H,niog(|H,|)r^, 

where R denotes an e-covering of A, and where the infimum is taken on all s- 
coverings. The (8, "f)-dimensional Hausdorff measure of A is 

(51) M^’"<{A) = lim M^’^{A). 


Since is composed of ~ randomly located intervals of length 3-2 , 

standard ubiquity arguments (such as in [T2[7]) yield that 


> 0 ; 


(491 implies that (which can be rewritten as a countable union) has a 

vanishing (rjuj, 2)-dimensional Hausdorff measure. Thus 


_ (J E^,\>0. 

\ 0J'<UJ / 

Since this set is included in we obtain that 

dim(il^) > Tjuj. 

It suffices now to rewrite these dimensions as a function of the p-exponent to 
obtain Theorem [TJ 


Numerical examples for the estimation of d^{H) and d^{H) of a lacunary wavelet 
series are given in Figure As predicted by theory, the numerical estimates of the 
p-exponent multifractal spectra are not invariant with p but follow the evolution 
with p of the theoretical spectra d'P{H). The positions of the mode of the estimated 
spectra have a constant negative bias; yet, quantitatively, they very well reproduce 
the shift of the mode of the theoretical spectra to smaller values of H for increasing 
p, revealing the lacunary nature of the function. A refined analysis is possible with 
the estimated lacunarity exponent multifractal spectrum d^{H), which has been 
computed here for several values of p for illustration purposes. The mode of the 
spectrum is estimated at H ^ 0.2 (instead of the theoretical H = 0.25). This 
clearly indicates the existence of positive lacunarity exponents. While the esti¬ 
mates for small values of p fall short of revealing the full support of the theoretical 
multifractal spectrum, they still enable one to identify a relatively large interval of 
positive lacunarity exponent values. The best estimate of d^{H) is obtained for the 
canonical value p = po = -boo {q = go = 0) in this example and produces a satisfac¬ 
tory concave envelope of the theoretical multifractal spectrum that provides clear 
evidence for ensembles of lacunary singularities with a range of positive exponents. 


5. Concluding remarks 

The analysis that we developed is based on the assumption that po(/) > 0, or 
that p/(p) > 0 for p small enough, so that p-exponents can be defined, at least, 
for p < Po; we saw that this assumption allows us to deal with distributions of 
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-©- p = I -■- p = 1.2 -- p = 1.7 -B- p = 2.5 -0- p = 5 - % -p = oo 


Figure 5. Lacunary wavelet series: A typical sample path of 
a lacunary wavelet series (a = 0.3, rj = 0.8, top row) and estimated 
structure functions (center row) and multifractal spectra (bottom 
row) for p-exponents (left column) and lacunarity exponents (right 
column) obtained with different values of p. The dashed lines in¬ 
dicate the theoretical multifractal spectra. 


arbitrarily large order and, equivalently, to model pointwise singularities with ar¬ 
bitrarily large negative exponent. However, this does not imply that any tempered 
distribution satisfies these assumptions. Simple counterexamples are supplied by 
the Gaussian fractional noises Ba for a < 0 whose sample paths can be seen as 
fractional derivatives of order ^ — a of the sample paths of a Brownian motion on 
K (Gaussian white noise corresponds to a = —1/2, in which case it is a derivative, 
in the sense of distributions, of Brownian motion). In |18| the wavelet and leader 
scaling functions are derived, and it is proved that rjB^ = —op, hence always is 
negative. However, the following result shows that, as soon as the wavelet expan¬ 
sion of the data has some sparsity, then this phenomenon no more occurs, and po is 
always strictly positive (note that this situation is quite common in practice since 
sparse wavelet expansions are often met in applications). 
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Definition 12. A wavelet series k *5 sparse if there exist C > 0 and 

t] < 1 such that, on any interval [1,1 + 1], 

Card{k : Cj-^ 7^ 0} < . 


Typical examples of sparse wavelet series are supplied by lacunary wavelet series 
or by the measures which satisfy (29). The following proposition implies that 
multifractal analysis based on p-exponents is always possible for data with a sparse 
wavelet expansion. 


Proposition 5.1. Let f be a tempered distribution, which has a sparse wavelet 
expansion, then r]f{p) > 0 for p small enough, so that po(f) > 0. 


Proof: Since / is a tempered distribution, it has a finite order, and thus it is a 
derivative of order A of a continuous function. Therefore / belongs to so 

that 


\cj,k\ < C2^^- 


Using again com pactly supported wavelets, the same argument as in the proof of 
Proposition 3.4 yields that there are at most 02"^^ nonvanishing wavelet coefficients 
at scale j; it follows that 


2-^ |cyfe|P < C2-^2^^2^P^ 
k 


so that rif{p) > 1 — p — Ap, and r]f{p) > 0 for p < (1 — ri)/A. 
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